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ABSTRACT 


The gas dynamic structures of the transport shock and the downstream 
collisional relaxation layer cire evciluated for partially ionized monatomic 
gases. Elastic and inelastic collisional nonequilibrium effects are taken 
into consideration. We account in our microscopic model of the atom for 
three electronic levels. Nonequilibrium processes with respect to popula- 
tion of levels and species plus temperature are considered herein. By 
using an asymptotic technique the shock morphology is found on a continuum 
flow basis . The asymptotic procedure gives two distinct layers in which 
the nonequilibrium effects to be considered are different. A transport 
shock appears as the inner solution to an outer collisional relaxation 
layer in which the gas reaches local equilibrium. The results show fo\or 
main interesting points: (i) on structiiring the transport shock, ioniza- 

tion and excitation rates must be included in the formulation, for the 
asymptotic method does not give frozen flow with respect to the population 
of the ionized and excited levels, (ii) the sharp rise in electron tempera- 
ture that might be expected after the transport shock is diffused to its 
beginning, (iii) the collisional relaxation layer is reduced rationally 
and accurately to qiiadrature for special initial conditions which (iv) are 
obtained from new Rankine-Hugoniot relations for the inner shock. A family 
of numerical examples is displayed for different flow regimes. Argon and 
helium models are used in these examples. 
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I. INTRODUCTION 

The possibility of obtaining large gradients by compressing a gcis 
through a shock wave has led many researchers to study transport, transfer 
and relaxation phenomena both theoretically and experimentally. 

It is well known that as the strength of the shock wave becomes lcu?ger 
and larger, radiative and collisional nonequilibrium effects begin to domi- 
nate and the thickness of the shock is determined by the characteristic 

lengths arising from such nonequilibrium processes. 

1-3 

Severed authors have devoted their attention to the relaxation 
zone behind a shock wave, taking into account different nonequilibrium 
phenomena to study the ionization behind it. An excellent review in this 
field is provided in Ref. 4. Jaffrin^ studied the transport shock in 
argon considering the degree of ionization frozen and without including 
an excited level in the atom model. Chubb^ included some transport 
properties and used a bimodal Mott-Smith velocity distribution function 
to structure the transport shock and the relaxation zone behind it. 

7 

Clarke and Ferrari studied the nonequilibrium ionization due to 

radiation and collision; they found a photoionized prec\u?sor in the shock 

0 

structure. More recently, Clarke eind Onorato structured shock waves in 
monatomic gases using asymptotic arguments. They found different asympto- 
tically embedded layers in which the nonequilibrium effects due to collisional 
and radiative ionization are uncoupled in each layer. Further applications 
extending this method are found in works by Pirri and Clarke,^ and Fainsworth 
and Clarke . 

Finailly, Foley and Clarke^^ refined the micTOScopic gas model of Ref. 8. 
They found that the shock morphology consisted of certain asymptotically 



2 


embedded layers: (i) a nonequilibrium precursor characterized by the 
absorption of radiation, (ii) a chemiccQ.ly frozen and transparent transport 
discontinuity, (iii) an optically transparent collisional layer characte- 
rized by elastic and inelastic collisions driving the gas to equilibrium, 
and (iv) a hot radiating tail in local equilibrium. 

The purpose of this work is to study the transport shock and the 
collisional layer in any partially ionized and excited monatomic gas. It 
will be shown, using an asymptotic argument, that although both layers can 
be separated asymptotically, they must be evaluated together, for they are 
intimately interrelated. Furthermore our results will show that the 
transport shock is not chemically frozen as considered by other authors: 
Therefore ionization and excitation rates need to be included in the 
formulation. 

Numerical examples are performed for models of argon and helium by 
using a continuum or Navier -Stokes approach, for different flow regimes. 
Thermal and chemical nonequilibrium effects are included in the formulation. 
The results of these calculations justify, a posteriori , the ideas 
introduced in the present work. 
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II. MODEL AND STATISTICS 

We consider the steady , one dimensioncil flow of a partially ionized 
and excited monatomic gas. No applied external electric or magnetic 
fields are considered herein. 

A. Microscopic Model 

The monatomic gas atom A is modeled microscopically by three 
different electronic levels, i.e. a ground level 1 which will provide 
species A^ ^ , an excited level providing A^^^ , and the continuum 
^ giving the single ionized ion Therefore, the gas is considered 

to be a mixture of atoms (both in the ground and excited levels), ions, 
and electrons. The term heavy particles, denoted by subscript h , is used 
when we refer collectively to atoms and ions. 

Since the mass of atans and ions is much greater than that of the 
electrons, the elastic energy transfer between heavy and light species is 
considered very difficult and therefore highly inefficient. 

We also consider only gases that are difficult to ionize, i.e. the 
ionization energy is larger than the translational energy throughout the 
entire flow process. The creation of an ion by inelastic collisions will 
result in the loss of an amount of energy equal to the ionization energy 
by the heavy gas or the electron gas depending on which species collides 
with the atom in the act of ionization. Considering the electron gas, we 
have two difficult processes which oppose one another: (i) Through elastic 
collisions the heavy particles and the electrons will try to equilibrate 
their translational energy, and (ii) every act of ionization by electron 
impact will substantially reduce the energy of the electrons. 
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Our gas model also provides for an excited level which is assumed to 
be very close to the continuum. Therefore excitation also requires a con- 
siderable amount of energy but ionization from the excited level is then 
possible with a much smaller amount of energy. 

As the gas enters the shock, the transport mechanisms will convert the 
initial kinetic energy into thermal energy. However, since the elastic 
energy transfer between heavy particles and electrons is very inefficient , 
heavy and light species will have different translational energy. A subse- 
quent collisional reloucation zone will then appear in which all species come 
to equilibrium. 

B, Statistics 

Elastic and inelastic collisions are considered in the present work. 

Mathematical expressions are then needed for the cross sections for the 

different types of collisions. Although few data are available for excited 

atoms, the statistics used generally are not the best we could hope for. 

There are two types of elastic collisions: (i) elastic collisions 

( 2 ) 

between species M and A , and (ii) elastic collisions between M and 

A^^^ , M denoting any species A^^^ , A^^^ , A^^^ , or electrons. It 

is known that the effective diameter for the former type of collisions is 

larger than that for the latter type of collisions. However, since few 

particles are in the excited level, the excitation energy being almost as 

large as the ionization energy, the effect of the additional size of the 

excited atom cross section should not be appreciable in the transport 
12 13 

properties. ’ Therefore, it will be assumed that the cross sections for 

( 2 ) 

elastic collisions between M and A are equal to the cross sections 
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for collisions between M and . 

More specifically, the elastic collisions considered are: collisions 
between heavy particles , collisions between electa?ons , and collisions be- 
tween electrons and heavy species. Since collisions between particles of 
the same species take place at a much faster rate than collisions between 
electrons and heavy particles, and the energy transfer between these two 
species is very inefficient, it is possible to define a heavy particle 
temperature Tj^ and an electron temperature in the thermodynamic 

sense throughout the shock. Therefore, will be, in general , different 

from T 

e 

For the inelastic collisions, we only consider transitions by electron 
impact. This is due to the fact that, since the gas is already sufficiently 
ionized ahead of the shock , the most efficient mechanisms for further exci- 
tation or ionization are those by electix>n-atom inelastic collisions. It 
has been shown previously that the treuisport shock and the collisional 
relaxation layer are transparent to any radiation. Therefore any radia- 
tive effects need not be included in the present work. Thus, the following 
parallel chemical reactions will account for the inelastic transitions: 


a'^> 

A«) 


+ e ^ A 
t e j A 
+ e ^ A 


(+) 

( 2 ) 

(+) 


+ e + e 
+ e 

+ e + e 


(1-c transitions) 
(1-2 transitions) 
(2-c transitions) 


( 2 . 1 ) 


As a result of the aforementioned nonequilibrium rate processes , the gas 
is not generally in local equilibrium with respect to population of levels 
and species or temperature. 

To make the problem more tractable mathematically, we now introduce 
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several £ priori assumptions which will be justified a posteriori : 

(a) Small ion-slip: The ion-atom elastic collision cross section due 

to the charge exchange mechanism is very large , so that the ion-atom mean 
free path is small compared to the shock thickness. There will then be 
many collisions between atoms and ions and therefore their temperature and 

5 

velocities will be equal. 

(b) Quasi-charge neutrality: Charge separation effects have been 

neglected because the Debye length is much smaller than any other charac- 
teristic length.^ Also the induced electric field is negligible.^ Thus, 
we assume that the ion and electron number densities are equal and there- 
fore all the species have the same macroscopic velocity. 

(c) Initially frozen elastic and inelastic collisions; It is 
possible for the upstream flow ahead of the shock to be locally out of 
equilibrium with respect to population of levels and species and/or 
temperature. In that situation a precursor will appear. This nonequili- 
brium precursor can be structured, for instance, by the absorption of 
radiation coming from a hot layer behind the shock. Since the characteris- 
tic length of this precursor is much larger than any characteristic length 
of interest to us,^^ the upstream flow is seen as frozen on the length 

sccile of the shock. Even in the case when T, and T are different 

h e 

ahead of the shock, small compared to its value behind the 

transport shock. Since the elastic collision energy transfer is propor- 
tional to that difference, we can create an artificial state of equilibrium 
upstream of the shock without introducing an appreciable error. This is 
done by subtracting from the elastic collision energy transfer term, 
denoted by E , its numerical value at the front of the shock E(a) . 

The inelastic collisions can be treated in the same way, since the electron 
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temperature upstream is sufficiently low to neglect the ionization and 
excitation by electron impact ahead of the shock. This assumption 
will be clcirified later on. 
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III. FORMULATION 

Under the assumptions and the model as described in Section II, the 
equations governing our problem are then the equations of state , continuity , 
momentum, energy, electron energy, ionization and excitation rates. Let 
n^^ , n 2 , n^ denote the number densities of atoms in the ground level, 
atoms in the excited level, and electrons (or ions) respectively, u the 
velocity, p the pressure, p the mass density of the mixtiare, a the 
degree of ionization, 6 the degree of excitation, the ionization 

temperature, the excitation temperature, R the gas constant per 

unit mass of atoms, and m^ the mass of an atom (or ion). Then we have 


p = m, (n,+n_+n ) , o = 
n X z e 


n 


n. 




6 = 


n, +n.+n 
1 2 e 


(3.1) 


The enthalpy per unit mass of the mixture, H , shall be written as 


H = R 


^T,+aT^) + aT 
2 h e 




(3.2) 


Then our equations are 


p = pR(Tj^taT^) , 


(3.3) 


pu = G , 


(3.4) 


. 2 4 du 

P + - 3 ’ 


(3.5) 


u^ . 5. 


*** ' ' u X 

+ ^(T +oT ) + aRT. + 6RT, - ii-K 

? 2 h e ] £ 3 dx ^ dx e dx 


= G 


W 

T ’ 


(3.6) 
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wherein (3.5) and (3.6) u , are the coefficients of shear vis- 

«'osity of the mixture, thermal conductivity of the heavy particles, and 
thermal conductivity of the electrons respectively. The electron energy 
equation can be written as 

3 ^^e ^ ^ _ 1 d j-w (3.7) 

2 dx u dx GRa dx ^ e dx ■' a dx a dx ’ 

E being the energy transfer term between heavy particles and electrons by 
elastic collisions, and the othtsr two relevant terms on the right denote 
the energy transfer between the same two species by inelastic collisions , 
leading to essential new effects . 

Let (k 1 . . denote the reaction rate coefficients in the chemical 
'■ mn-' 1- j 

reaction between particles of the mth and nth species for transition from 
the ith to the jth level. Since as discussed in Section II, the only 
reactions considered are those by electron impact , we let 


Tk 1 • . — i A ~~ A 

mn-'i-j ^ eh'' 1-3 1-3 


(3.8) 


Then from (2.1) and (3.1), the population rate equations are 


da „-l 

^ "'v, 

dx h 




(3.9) 


^ _ ..-1 

dx 


= G m^ 




(3.10) 

jf inelastic collisions by heavy— heavy encounters were important in 


r 
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the ionization and excitation processes , more terms would be required on 
the right hand side of the rate equations. According to our assumptions, 
transitions by heavy particle impact are negligible and therefore are not 
included in the formulation. 

It is expeditious to utilize the fact that the reverse reaction rate 
coefficients ic^ ^ can be expressed as functions of the different forward 
reaction rate coefficients ^ through equilibrium considerations. Then, 
from (3.9) and (3.10) we obtain 


c-1 e 1-c e p 2 


'^c-2^’’e^ ■ '^2-c^'^e^ 7“ ~ ’ 

“e 


(3.11) 


’'2-1^'^e^ f '^i- 2 ^’’e^ 


^■“e'^e 


(3.12) 


where the subscript E denotes reference equilibrium conditions at p 

and . The remaining k's can conveniently be written as lengths for 

7 

gas dynamic purposes. A local characteristic length for the ionization 
process is seen to be 


X 

a 


V 

p(l-a-B) 


'1-c 


(3.13) 


Similarly, from Ref. 5, we have expressions for the different transport 
coefficients and the elastic collision energy transfer term E expressed 
in terms of lengths 

y = p(RTj^)^^^X^ , (3.14) 

= pR(RTj^)^^^Xjg^ , (3.15) 

= apR(RT^)^^^Xj^^ , (3.16) 
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E = (RT 1 ^ (3.17) 

6 U A 

e 

where , Xj^^ , Xj^^ , and X^ are lengths characterizing the different 
tramsport and transfer phenomena. In the Appendix expressions are provided 
for the different forward reaction rate coefficients and for the character- 
istic lengths defined aibove as functions of the dependent variables. Thus 

we have seven equations in the unknowns u, p,p,T, ,T ,a, and 0 . 

h e 

The equations will now be written in dimensionless form. We define 
a dimensionless independent variable n through the relation 


X 

dn = X ^dx , n = X ^(x)dx . (3.18) 

a I a 

■’o 

Let be the electron temperature at station b inside the shock. 

The significance of station b will become clear later on. Then the 
following dimensionless dependent variables are found natural and appro- 
priate to the problem 


T. 


U 


u 


e'b 


CR(T^), ] 
e h 


172 


[R(T ) 
e b 

, p = p , p 


G[R(T 


-^ 1 7 2 .( 3 - 19 ) 
e>b^ ^ 


where T^ stands for T^ , T^ , T^ , or T^ accordingly. 

We substitute the relevant lengths for the transport coefficients and 
the elastic collision energy transfer term into the equations and make them 
dimensionless following (3.18), (3.19). We have 

A A A A 

p = p(Tj^+oT^) , (3.20) 


pu = 1 




(3.21) 



NJ| cn 


12 


; + u i ^ ^ _ V 

P - 3 x„ o\ dn - „„ ) ,1/2 • 

e b 


(3.22) 


(VaTe) + aT . + 3T^ + ^ T, 


^ 4 ;i/2 du ^Kh -1/2 ‘^’’h 


dn 


pT 


dn 


^Ke —1/2 

X dn " 2R(T ) ’ 

a e D 


(3.23) 


o dT 
3_ e 

2 dn 

3 ^ 


X 

a dn ^X 


-1 ^ I A fDA 

Cl dn ~ « 'I" 


a 


dT 
e 

dn 


X t^/2 

-i ) 

3 


— ^ 
o dn 


T 

1 Z dB 


a dn ’ 

3 


(3.24) 


i. ^ 

a dn 


2 l-a„-B_ 

•Co ^ 

2 Bt. 

, a E 

1 2-c a 

1 - “ 

1-0-6 2 

K, _ 1-a-B 


L “E J 

1-c — 

a 

L “d 


= I 


1* 






1 dB _ ''1-2 

. B 

1^2 c ^ 


a dn K, ^ 
x-c 

1-a-B Bj. 

"im: 

a 

E 


= I. 


(3.25) 


(3.26) 


From the above equations, we see that we have four ratios of character- 
istic lengths and two ratios of forward rate coefficients which could be 
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converted in length ratios very easily, but we do not wish to do so here. 

Now, from the definition of X and X^ , we can see that both 

e Ke 

lengths are comparable and of the same order of magnitude. It is also 

clear from the expressions of X and X,,. that both are of the same 

y Kh 

order of magnitude. Moreover, numerical considerations tell us that X„ 

Ke 

and Xj^j^ can be compcirable and that X^ is the largest of all the lengths 
introduced herein. Usually, the above lengths satisfy the following 
inequalities : 


^y ^Kh ^ ^Ke 


= X « X 
e o 


Notice that since the excited and the continuum levels are considered 
to be very close together, X^ also characteriz® the excitation process. 

From the above ordering of the characteristic lengths we see that if 
the problem or phenomena of interest to us occurs in a length of order 
X^ , we can neglect all the transport effects in (3.20)-(3.26). 
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IV. ASYMPTOTICS 


Let us now construct a mathematical repi^esentation of our physical 

model discussed in Section II-A. The gas can be represented by constructing 

thi^ee limits. If m is the electron mass, let 

e 



I 



(4.1a ,b,c) 


such that 




= 0(1) , s i 5 ^- = 0(1) , Fj s 


_ 2-c 


T. 




= 0 ( 1 ) , 


(4.2a,b ,c) 


the reeisons behind (4.2) will appear clear as we proceed. 

Notice that only one of the three limits (4.1) introduced is indepen- 

A 

dent. The last two are related through (4.2a), and 6 and T^ are 
connected in . In F^ an additional relation is provided between two 
of the three reaction rate coefficients. 

Ft*om the physics of Section II, and from (4.1c) it follows that 


f - 0 , (4.3) 

which means that the excited level is highly depopulated, i.e. because the 
excited and the continuum levels are very close together, it is easier to 
ionize from the excited level, rather than from the ground level, resulting 
in its consequent depletion. Similarly, from our physical model and (4.1c) 


it follows that 
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1-2 2-r 

-j^=0(l), (4.4a,b) 

K. K, 

1-c 1-c 

From our expressions for the forward reaction rate coefficients (Appen- 
dix B) and our limit construct (4.1) together with (4.2a) we can verify 
the validity of (4.4). Notice that (4.3) and (4.4b) are not independent 
for they are related by (4.2c). 

At this point, under our limit, other dimensionless dependent 
variables, e.g. , u , p , and p are considered 0(1). We are 

now in a position to study our (5quations and to try to simplify them by 
using the adaove orders of magnitude under the limits constructed. 

A. Outer Solution 

First consider the collisional relaxation process. The characteristic 

length scaling such phenomena is seen to be X defined in (3.13). We 

a 

assume that derivatives with respect to n ''' x/X are 0(1) under our limit. 

a 

From the discussion of the different characteristic lengths introduced in 
Section III and from (4.2b), we see that X^ is much larger than any of 
the other characteristic lengths. Therefore, on studying the collisional 
relaxation, it is permissible to neglect all the transport effects in (3.20)- 
(3.26). Thus, the layer in which the transport mechanisms are important, 
referred to as the transport shock, is seen to have zero thickness on the 
length scale defined by X^ . The characteristic electron temperature 
(T^)j^ introduced in Section III appears now as the electron temperature 
just behind the transport shock. Therefore in the outer solution or colli- 
sional relaxation layer, the momentum and energy equations (3.22) and (3.23) 
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become algebraic! It is possible then to evaluate the integration constants 
in the right members of (3.22), (3.23) anywhere inside the collisional layer. 
Evaluating them just behind the transport shock at stage b, (3.23) becomes 


2^'^h'’‘“'^e“^\^b““b^ + T.(o-aj^) 


1 + 


T. “-“b 


^2 *^2 
u -u 

+ 0 


(4.5) 


From (4.1c) and (4.3) we can neglect the second term in brackets in 
(4.5) as compared to one. Since the first and last terms are 0(1) , it 
is necessary that 


Tj(a-ab) = 0(1) » or “-Ojj = o(-t^). (4.6) 

*^j 

Our limit (4.1b) strictly implies that a does not change, from its 

value behind the transport shock, over the collisional layer. However, it 
• 7 8 

is known ’ that, for appreciable ionization, the electron temperature does 
not have to be close to . For all the numerical examples considered , 

^ A 

tends very weakly to its limits (T^.-vlO) , and therefore considerable 
ionization can occur. Moreover, by neglecting (1/T^) compared to unity, 
the ionization can be predicted very accurately with a 16-20% error. 

We now turn to the population rate equations. The brackets in (3.25), 
(3.26) account for the combined forward and reverse reactions. Their 
numerical values vary between zero and one. Only when the flow reaches 
equilibrium, will the forward and reverse reactions balance and give zero 
for the bracket terms. 

In (3.25). both the right and left members are 0(1) under our limit. 
Then, the population rate equation for the ionized level remains differential 
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over the outer solution. 

In (3.26), the left member is 0(p/a) which is small under (4.3). 

The right hand side is 0(1) under (4.2c) and (4.4a). Thus, the left 
member of (3.26) is negligible and 8 is given locally in terms of the 
other dependent variables by equating to zero the right-hand side of 
(3.26). The population rate equation for the excited level reduces to 
algebraic! 

If inu(3.24) we replace (!5ee Appendix) and by their respec- 

A 

tive expressions and , we see that the right member is 0(1) 

while the left-hand side is and therefore it can be neglected 

T. 

1 

under (4.1b). The electron energy equation also reduces to algebraic over 
the outer solution. Namely, it reads 

E - = 0 , (4.7) 

which means, physically, that the heating (cooling) of the electron gas 
through elastic collisions with the heavy particles is balanced by the 
cooling (heating) due to inelastic collisions between the same two species. 

We have shown that, over the outer solution, the system of five 
differential equations (3.20) - (3.26) reduces to only one differential 
equation, the one governing the population of the ionized level. Then, all 
the dependent variables, except the degree of ionization, will adjust to 
local conditions and can be expressed as functions of a , the only non- 
equilibrium variable, since a^, and 8^, can be expressed in terms of 
T^ and p . We say that, over the collisional layer, the flow is in a 
state of "autolocal equilibrium." The outer solution is then very simple 
to obtain, for the problem is reduced rationally to quadrature. 
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B . Inner Solution 


By reducing our system of equations over the outer solution to only 
one differential equation , it is clear that the boundary conditions at the 
front of the shock wave, station a, cannot be satisfied. On the scale X 

a ’ 

this means that a geometrically thin boundary layer must exist in which at 
least some of the derivatives in (3. 22)- (3. 24) and (3.26) tend to infinity. 
Considering the two lengths appearing in (4.2b), x must be scaled by X 

e 

instead of X^ according to 


dc = X ^dx 
e 


or 


? = 


X 

* 

X"^dx . 
e 
o 


(4.8) 


Then over the boundary layer we obtain 


A A 


p = P(V«T^) , 


pu = 1 , 


(4.9) 

(4.10) 


n + n - i _1L du _ V 

^ ^ [R(T ) ]^^^ ’ 

e b 


(4.11) 


|(Tj,+aT^) + aT^ + 6T, 


e 



4 \ - 1/2 

3 X h dc 
e ^ 


^Kh ;;i/2 ‘^^h 

— dT 

e ^ 


^Ke -1/2 ^e _ 

X^ ^ e d^ " 2R(T ), » 

e e b 


(4.12) 
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3 ^ 4 . ^ _ 1 A 

2 ^ dc ' a dc 


^Ke ;-l/2 

— dT 


ml/ 2 

T 

e 

u 


(T. -T ) 
h e 


_ i. [t — + T 

a Qj dc ^ dy ’ 


1 ^ 
a d^ 


A A 

T. ^1 
1 


1 - 


1-a-e 







+ r 



2 

* ^2 

^ " e 2 


“e -j 


“eJ 

4 


(4.13) 


(4.14) 


1 ds _ 1 1 

•^1-2 

3 ^'V^E 

- r 

Bp 

1 - “ ^ 


« d? m I’ 1 

^ T. 1 

3 

•'l-c 

k 

l-a-3 3g 

^2 

E 



(4.15) 


It is clear that in a layer scaled by « X^ , no appreciable 

ionization or excitation can occur, since X characterizes both the 

a 

ionization and excitation processes. Effectively (4.14) states that 


— ^ is which is small under (4.1b); however is o(l) . 

a d? ad? 

j 

The right-hand side of (4.15) is o(-^) while its left member is higher or- 

T. 

3 

der because everywhere 3 is' much smaller than o . Therefore in (4.13) we 

could neglect “ ^ as compared to “ ^ • f’or the reasons explained in 
the Appendix we prefer to substitute these derivatives by their corresponding 
expressions in the right of (4.14) and (4.15). Since X,, and X are 
comparable lengths, the left member of (4.13) is 0(1) . Therefore, although 
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ionization and excitation rates are small, these rates must be included in 
(4.13). From the ntimerical point of view, it is possible to consider a 
and 6 as constant at their Initial veilue over the inner solution; therefore, 
we do not wish to carry the asymptotic argument any further. 

Notice that, since can be con^arable to , Xj^^ and X^ the 

transport terras must be included in the momentum and energy equations. In 

A A 

(4.12), oT. and 3T will csmcel with similar terms appearing in the 
integration constant W . 

Then over the inner solution, (i) the transport effects are relevant, 
and (ii) although ionization and excitation are small, the inner solution 
does not give frozen flow with respect to the population of ionized and 
excited levels, for excitation and ionization rates must be included in the 
electron energy equation. We shall take a and g constant at their 
initial value. We say that inside the transport shock the flow is in a 
state of "slight ionization and excitation." 


Finally, we must mention that the error factor introduced by the 


asymptotic procedure is seen to be 



over both the outer and 


inner solutions. 
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V. JUMP CONDITIONS 

It is possible for the upstream flow not to be in local equilibrium. 

In that case a precursor will appear ahead of the shock front (see Section 
II). However, the flow in this prectirsor is seen as frozen in the length 
scales defined by either or . There is also a possibility, 

depending on the flow conditions , for the existence of a hot radiating 
tail in equilibrium, behind the colli si onal layer, in which the flow will 
reach complete thermodynamic equilibrium.^^ Even in that situation, since 
this tail is structvired by the absorption of radiation, the flow in this 
tail is seen to be in complete equilibrium on the length scale defined by 
Aq . It is possible then, from our equations to evaluate the jump condi- 
tions, denoted by [ ] , across both the outer and inner solutions. 

A, Outer Shock Jump Conditions 

At the end of the collisioncd layer (station c ) , thermal equilibrium 
between the species is reached, so that T, = T = T . Then, from our 
equations (3.20) - (3.26) we can get the outer jump conditions, because the 
flow derivatives outside the shock vanish in the length scale of the outer 
solution. Therefore we get, in the physiceuL variables 

a 

= 0 , (5.1) 

c 

together with the equation of state and 


pu 


= 0 


[p . puf . 


2 


(T, ) = (T ) = T 

he e c c 


a = 
c 




«c = 6 e«c’'’=> 


(5.2) 


Therefore, the flow at the end of the collisional layer will be in local 
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equilibrium. 


B. Jump Conditions Across Inner Shock 


Since « X^ , the flow derivatives in the collisional layer 
measured in the length scale will be vanishing small. The same state- 

ment holds if there is a precux>sor ahead of the shock. Then, if we know 
the flow conditions ahead of the shock (station a^) we can evaluate the 
flow conditions just behind the transport shock (station b) from (4.9) - 
(4.15). These jump conditions cire given, in the physical variables, by 



together with a=a^=Oj^ , equation of state and 

E(b) = [X^(b)r^[T^yb)+T^l2(b)] , (5.4) 


from the electron equation (4.13). 

Note that (5.4) does not give the generally accepted 



0 across 


the transport shock, but tells that there is an energy balance between the 
electron-heavy particles elastic and inelastic collisions. 

Although the conditions of shock discussed here do not appear to be 
classical Rankine-Hugoniot relations, they effectively are. For the type 
of shocks considered and under the assumptions made, (5.4) plays the role of 
a thermo<fynamic equation of state. 
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VI. METHOD OF SOLUTION 

From our previous discussion, the method for structuring the shock 
appears clear: (a) Evaluate the jump conditions across the inner and outer 
shocks as discussed in Section V. (b) Solve the collisional layer: the 

^ 1C I *7 

forward reaction rate coefficients can be expressed as functions of ’ 

and Qj, and are also functions of p and T^ . Then from (3.20) - 

(3.24) and (3.26), which reduce to algebraic form over the outer solution, 

A A A A ^ 

we can express p, u , p , T^ , T^ and 6 as functions of a alone. 
Substituting these functions in (3.25), we have an equation of the form 

= <'''• 

Eq. (6.1) is integrated numerically, by Simpson's rule, between the 
limits a=a^=Oj^ and which are known from the jump conditions. This 

integration will give a=a(n)=o(x) . Once a(x) is known, we can con 5 >ute 
the spatial distribution of the remaining dependent variables since they 
are given in terms of a . The solution over the collisional layer is very 
simple to obtadn because the problem is reduced to quadrature, (c) Solve 
the inner solution: The constants of integration are evaluated at a . 

Considering a and 6 constant over the transport shock, we introduce an 
auxiliai^ dependent varicible z through 



u 


^Ke ^ 

e 


( 6 . 2 ) 


Using the equations of state and continuity to eliminate pressure and 
density, we solve for the different derivative terms. We then eliminate 
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the independent variable ? , by dividing all the equation by (6.2) to 
get 


T +aT (T, ) +o(T ) 
he ha e a 


+ u - u 


du 

dT 


u ^ T, 

^ V ( h-v 

^ K K ^ 

Ke T 


1/2 


(6.3) 


dT. 


dT 


I <W> 




1." " .2 

^u-u^) - az 


1/2 


Kh r h 




Ke T 


(6.4) 


dz 

dT 


3 ^e ^ X II ^ 

2 ^Ke ^1/2 " " ; dc 


jl/2 

e 


u 


( Tv - T „) 
h e 


- Y 


T. 

da 

a 

dC 


— f do- j 
dr " IdC-' 


dc 


f— 1 


-r-1 


u 


K ^1/2 

Ke T 

e 


(6.5) 


wherein 


"a = TTt 




e a u 


Ph>a-<^Q - 


( 6 . 6 ) 


and in (6.5) ^ ^ ^ should be replaced by their appropriate expres- 

sions (6.3), (4.14) and (4.15) in terms of the other variables. 
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Numerical integration of (6.3) - (6.5) is initiated at the four- 
dimensional singular point at state ^ using the limiting values of the 
du '^'^h dz 

derivatives — ^ as determined from (6.3) - (6.5) by 

dT dT dT 

e e e 

I'Hdpital's rule. 

The numerical integration scheme used is a modification of the method 
described in Ref. 18. Eqs. (6.3) - (6.5), in vector form, can be written 


^ = F(w,T ) , (6.7) 

dT ® 


where w = w(u,T^,z) . Then 


w 


n+1 


= w + A • 

[ci] - ^A ]] 

-1 

F + 

(^-)1 

n 

1 2 n J 

1 

n 2 

3T 

e n 


( 6 . 8 ) 


where A is the step size, tl] the identity matrix, and [A..] = 

1 ] n 

3fi'| 

, i,j=l,2,3. Matrix [A] is evaluated numerically at each step. 

'■aw.-* 

J n 

To obtain the spatial distribution of the variables , nvimerical integration 
of (6.2) is required. Once the spatiaJ. distribution is known, the induced 
electric field and charge separation can be computed.^ 

The computations , which were performed on the Brown University IBM 
360 model 67 computer, require less than 5 minutes to perform. 
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VII. RESULTS AND DISCUSSION 

Numerical examples have been calculated for microscopic models of 
argon and helium. The data are chosen from Ref. 5 (Figs. 1,2) and Ref. 11 
(Figs. 3-10). All the results obtciined for the collisional relaxation layer 
and the transport shock are plotted as functions of the physical distance 

X . 

Both the inner solution (a to ^) and the complete solution (a to c) 
discussed here are true shock waves in the sense of nonlinear sound waves. 
Although the interaction between electrons and heavy particles is weak, the 
electrons are constrained by electrical forces and are forced to shock. 

The transport shock and the collisional layer are interrelated, since the 
former arises as the inner solution to the latter. Since the flow over the 
inner solution is not totally frozen with respect to population of the 
excited and ionized levels, ionization and excitation rates must be included 
over the inner solution. Inelastic collisions are responsible for cooling 
the electrons as the flow approaches the collisional layer, and they thus 
prevent the light and heavy species from reaching thermal equilibrium on 
the small scale of the transport shock, as would otherwise occur through the 
elastic collision mechanism.^ The inclusion of the inelastic collisions is 
also responsible for the new Rankine-Hugoniot state relation replacing 

= 0 across the transport shock. Although Foley and Clarke did not 

obtain the correct value of T^ behind an unstructured transport shock, 

the electron temperature catches up to its correct value very quickly at 

the beginning of the collisional layer for their numerical examples. In the 

context of their assumptions, their treatment of T^ across the transport 

shock is self consistent. However, an incorrect initial value in T for 

e 
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the structure of the collisionauL layer could give substantial discrepancies 
in other dependent variables (cf. Fig. 3). 

Our results for the outer solution were anticipated, on numerical 
grounds, in Refs. 4 and 15, although our treatment of 3 over the 
collisional layer is different from Ref. 4. In our model, no assumption is 
made that every excited atom created will be immediately ionized and thus 
a more reeilistic treatment of the ionization process is provided. Figs. 1 
and 3 show typical solutions over the collisional layer for the conditions 
noted. The asymptotic result is compared to the exact solution obtained 
by solving simultaneously Eqs. (3.20) - (3.26) neglecting the transport 
properties . 

As the transport phenomena are included in structuring the inner 
solution , an electron precvu?sor evidently appears due to the electron 
thermal conductivity (cf. Figs. 2, 4-10). The sharp rise in electron 
temperature that might be expected behind the transport shock as a conse- 
quence of the jump condition 1 t~T^ = 0 » actually occurs at its beginning 


Lej 


due to K 


e 

Over the inner solution, the actual rise in the degree of ionization, 
for all the examples evaluated, is less than 6%. In Fig. 2, the inner 
solution was computed twice for two disparate values of 3 . We can see 
that the two results are almost identical. This is due to the fact that 
3 is everywhere much smaller than a , its numerical value not affecting 
the results. Thus, our results clearly veuLidate the treatment of o and 
3 as constant at their initial values over the transport shock. 

If we are not interested in the transport effects, it is possible to 
disregard the structure of the inner solution. However, with the use of our 
jump conditions , we can evaluate the correct flow variables behind the 
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transport shock. One could neglect the structure of the heavy particle 
transport properties (cf. Figs. 5, 8, 10) and treat this as a discontinuity 
representing an "innermost" solution within our inner solution, structured 
solely by electron conductivity; but this simplification will not always be 
possible, because matching difficulties will appear. This innermost 
solution can not be realistically embedded by an asymptotic argument on the 
length scale of the inner solution. For the numerical examples evaluated, 
the thickness of this innermost solution is about 25-30% of the total thick- 
ness of the inner solutions given. 

There are some cases in which the asymptotic method discussed here will 
not give satisfactory results particularly in the case of the lighter gases, 
depending on how strongly the gas parameters will tend to their respective 
limits (4.1) and (4.2). If the gas is lighter (6 diminishes), the elastic 
coupling between electrons and heavy particles increases. For helium, for 
instance, only one of the numerical examples considered (cf. Fig. 6) was 
successful, the one corresponding to the smaller degree of ionization of 
the upstream flow (a = 0.0412). As a increases (a 'v- 0.1), the electron 
temperature behind the transport shock increases, due to coupling effect, 
by a factor of 2 with respect to a similar argon example. At these high 
temperatures, and with the initial concentration of electrons corresponding 
to = 0.1 , and X^ can become comparable and our asymptotic pro- 

cedure will fail, since appreciable ionization will occur inside the trans- 
port shock. The inner solution for the helium case we have been able to 
evaluate is shown in Fig. 6. The outer solution for this example can be 
found in Ref. 15. 

All the assumptions upon which this model is based are completely 
justified £ posteriori . The values obtained for the induced electric field 
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and charge separation are indeed very small. The initially frozen elastic 
and inelastic collision assumption used to structure the shock is also 
justified; a difference of 2 - 3% in the results appears according as it 
is used (as shown explicitly in Eq. (6.5)) or not used. 

19 

Finally, we must mention the work by Magretova, Pashchenko and Razier 
who studied a strong shock in air. Taking into account the thermal con- 
ductivity of the electrons as the only transport mechanism, they found an 
electron precursor. The heating of the electron gas is said to produce 
appreciable ionization before the heavy gas is shocked. With our model, 
no such ionization appears in the electron precursor for any of the exam- 
ples we have evaluated. 

The method presented herein can be used, if desired, together with 
the work of Foley and Clarke^^ to obtain the full structure of a shock 
wave when radiative as well as collisional effects are to be taken into 
account, providing a rational description of the kinetics of the collisional 
and radiative ionization reactions at high temperatures. 
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APPENDIX: ALTERNATE ARGUMENT 

Over the outer solution it is clear that — 4^ is much smaller than 

a dn 

1 dot 

— ^ as seen from (3.25) and (3.26) under the asymptotic argument in 
Section IV. Therefore, in (3.24) we can neglect the former term as compared 


to the latter. If in (3.24) we substitute t>y its corresponding 

expression I^ and neglect the left hand side which is 0 (-4^) , and ^ 


T. 

D 


we get 


1 ^ 
dn 


;i/2 


- 


u 


2 1-a^-e^ 


2 6., 

. a E E 

^ P 


l-a-6 2 

^2 


“r 

_ ' _ 


“e 


= 0 


(Ap.l) 


Now, from (3.26), we can write 


PfI 

1 - — — 

_ ’'1-2 

B ^"“e~^E 

^62 


1-a-B 

“e 

1-c 

E 


0 © 


(Ap.2) 


Thus (Ap.l) can alternately be written as 

' 1/2 
u 


2 1-a^-B^ 


„ l-a„-B 

a E E 

1-2 

, B E '^E 

1-a-B 2 


1-a-B Bp 

“e 

1-c 

ia 


0 .(Ap.3) 


We can verify very easily that (Ap.3) is precisely the form (3.24) takes on 
(besides a factor of 



T. T. 

3 3 


multiplying the last term) when we substitute in it — 4^ and — 4^ by 

a dn a dn 
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their corresponding expressions 1^^ and • Ths difference between 
(Ap.l) and (Ap.3) is seen to be of a higher order. 

On physical groxmds, we prefer Eq. (Ap.3) to Eq. (Ap.l), because 
(Ap.3) states explicitly that, over the outer solution, the three difficult 
processes introduced in Section II, namely energy transfer between heavy 
and light species through elastic collisions, ionization from the ground 
level, and excitation from the ground level give the energy balance of the 
electrons . 

The above argument also holds in the discussion of the inner solution 
since ^ ^ can be neglected as compared to ^ ^ , because 3 is every- 
where much smaller than o according to the microscopic model of Section 
II. We also prefer, for the reasons explained above, to substitute in 
(4.13) ^ and ^ by their expressions in the right of (4.14) and (4.15) 

respectively, to get 


3 ^ 
2 dc 


dc a dc 1 e dC 
e 



1 . 


^ l-a^-8^ 

A 

T K 

^ 1 1-2 

, 6 

I 

1 


1-a-B 2 

E _ 

1 K- 

T. 1-c 
3 

l-a-6 





( Ap . 4 ) 


It is clear that the last term in the right of (Ap.4) can be neglected 


under (4.1c). 
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Thus we see that in the outer solution as well as in the inner 
solution, substituting the ionization and excitation rates by their cor- 
responding expressions is equivauLent to neglect the excitation rate as 
compared to the ionization rate, the difference between these two altema 
tives being higher order. However, we prefer the former to the latter 
because of its physical mecining, particularly over the collisional layer. 


I 
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APPENDIX A 
TRANSPORT PROPERTIES 


The mathematical expressions for the different coefficients of 
viscosity and thermal conductivity are taken from Ref. 5. If denotes 

the elastic collision cross section for an encounter between particles of 
the ith and jth species and subscripts ^ ^ » ^^id ^ refer to atoms (in 
the ground and excited levels ) , electrons and ions respectively , we have 

Heavy transport coefficients 

The coefficient of shear viscosity of the mixture is 

, (A.l) 


where is a characteristic length defining the viscosity of the heavy 

particles. We can express in terms of the different elastic cross 

sections and therefore in terms of the dependent variables, i.e. 


X 

U 


_ Si/iT 1 
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a aa 


1 -t 
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e la 
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aa tt 
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(A. 2) 


or 
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~ 6 p 
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aa 
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1 + _2_ -il 
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^ 1-a a. 
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(A. 3) 


The coefficient of thermal conductivity of the heavy particles , is 
given by 


1 
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(A. 4) 


being a characteristic length for the heavy thermal conductivity and 
it is defined by 


Kh 


= 75/if 1 

" 64 p ^l-a)o 


aa 


1-a 


aa 


a. a. o. . 

1-a o 1-a o. 

aa la 


(A. 5) 


From the definitions of X and X„. 

U Kh 


we can see that 


^Kh _ 75 1 - 

X " 64 5 " 32 • 

y 


(A. 6) 


Thus both lengths characterizing the transport properties of the heavy 
gas are almost the same. 


Electron transport properties 

The viscosity of the electrons can be neglected as compared to the 
viscosity of the heavies due to the small mass ratio of the two species. 
Since y^ = 0 , the viscosity of the mixture will then be given by 

" = "h * “e = "h • <*•’> 


The electron thermal conductivity is given by 


(A. 8) 


where X^ is a characteristic length for the electron thermal conductivity, 
defined by 
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(A. 9) 


The energy transfer between heavy and light species through elastic 
collisions E is given by^ 


E = ^ , (A, 10) 


Xg being a characteristic length defining this type of energy transfer, 
defined by 


X = i g 

e 8 V 2 p ao . 

ei 


1 + 


I V/ 

-a ea 

a a . 
ei 


(A. 11) 


It is clear from (A. 9) and (A. 11) that both and X are 

Ke e 

comparable lengths. 

Moreover, on numericcil grounds, and depending mainly on the flow con- 
ditions, the chau?acteristic lengths defining the transport properties X 
Xj^ , and Xj^^ can be comparable, and therefore » ^Kh^\e 

non-neglibible quantities. Our results prove the Veilidity of the above 
argument . 

Since 


X 
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^Ke ^e 


X,. X 
Ke e 


= o(-;^) , 


T. 


(A. 12) 
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^Kh ^Ke K 

X X 
Ke e a 


= o{^) 


T. 

3 


(Aol3) 




(A. 14) 


We can neglect all the transport properties over the outer solution charac- 
terized by X^ within the orders of magnitude and limits constructed in 


(4.1) and (4.2). 
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APPENDIX B 

COLLISIONAL STATISTICS 


In this appendix expressions for all the cross sections and reaction 
rates are provided for both argon and helium. 


Elastic Cross Sections: 

1. Argon: All the cross sections for cU?gon are taken from Ref. 5. 

These cross sections in terms of the flow veiriables are: 
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(B.l) 


(B.3) 


where e is the magnitude of the electronic charge , and k the Boltzmann 
constant. The remaining cross sections are 
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C-0.35 + (0.775xl0"^)T ] x 10"^®cm^ 

e 


for T >10 ®K 
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y (B.4) 


[0.39 - (0.551x10"*^)! + (0. 595xlo"®)T^] x lo'^^cm^ for 

^ e e 

T < 10*^ °K 
e - 
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2. Helium: The cross sections for collisions between charged 

particles are the same as for argon, i.e., and a are given by 

(B.l) and (B.2), since the cross sections are independent of the masses 
of the particles involved. 

From Ref. 11 

a = 5 X lO'^^cm^ . (B.6) 

ea 


20 

From Hirchsfelder et al, the average cross section for the computa- 
tion of viscosity and thermal conductivity is given by 




f 


3 f 9 

X q' ■'(xkT)exp(-x)dx , 


(B.7) 


.( 2 ) 


where Q denotes the interaction potential for the colliding molecules. 

From Ref. 12, the interaction potential for an atom and an ion is 
given by the empirical expression: 


= [48.0-6.96iln(xkT+0.4) xio'^^cm^ (B.8) 

expressing kT in eV . Expanding the logarithmic terra 

iln(xkT.0.4) = l„0.4 4 2|ii^ 4 | ^ * — ] 

= anO.4 + 2 — . 


Substituting (B.9) into (B.7), we obtain 
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Since 


0 < 


T 



< 1 


we take an average for the above integral to get 


a . 

la 


71 X io"^^cm^ . 


(B.ll) 


From Monchick, the average cross section for atom-atom collisions 
is given by 
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8o p I 
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— * = 2a p I 
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(B.12) 


where a = An and A and p are related in the expression for the 
interaction potential 


= A exp(-r/p) , (B.13) 

with p in angstroms. Monchick gives a tabulation for I/,, as function 

of a . From Ref, 22 we have for the H -H interaction 

e e 

A = 157 ev, p = 1/3.93 A . (B.14) 


For our range of temperature 300° < \ < 130000®, we can extrapolate 
Monchick's results. The difference in the values of is not very 

significant however. By taking an intermediate value 

^(2,2) ~ 0*512267 , (B.15) 

corresponding to a = 4.5 , we obtain 

= 0.066291 [14.4137- An x io"^®cm^ . 
aa n 


(B.16) 
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Inelastic Processes 

All the inelastic processes considered are those involving chemical 
reactions by electron impact as given by (2.1). 

According to these reactions, we can write 


*^^^e 3 3 

= (C, n,n - K ,n + tc_ n.n - ic _n , 
dt 1-c 1 e c-1 e 2-c 2 e c-2 e ’ 
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where in (B.i?) and (B.18) the number densities can be written in terms of 
p , a cuid 8 . We know from equilibrium considerations that 


n. 


‘ 2-1 




= IC 


1-2 8 ^ ’ 


(B.19) 


= K, 


/■ l^ 


'c-1 1-C 


n 


= K 


"h 


1-C p 


(B.20) 


'c-2 ^2-c 




n 


= K 


2-C p 2 ’ 


(B.21) 


where E denotes equilibrium conditions. 


We wtite a^, and 6^ in terms of the dependent variables. Then 
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(B.24) 


where k is the Boltzmann constant, h the Planck constant, and g^ and 
gt are the degeneracy factors in the electronic partition functions for 
atoms and ions respectively. 

We now need expressions for the forward reaction rate coefficients. 
From Refs. 16 and 17 we have 
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where x is the potential difference between levels p and q 
p-q 

expressed in eV , and 


A = 




A+30 

^ ■ lOA+25 • 


(B.26) 


The numerical values taken for argon and helium are 


Helium 



^ = 285000°K 



T = 244800OK 
£ 


0 , 


(B.27) 

(B.28) 


and 
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gQ = 1 , = 12 , go = 2 

Argon 

T. = 183000°K , T. = 13402QOK , T = 2065®K . (B.30) 

(.3 3 P 

Expressions (B.22) - (B.24) are to be substituted in the rate equations 
for the population of the excited and ionized levels in the appropriate 
dimensionless form defined in Section III. 
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LIST OF CAPTIONS 


FIG. 1. Outer solution. Argon, = 1.5xl0~^ atm., (Tj^)^ = ('Tg)g = 1000°K , 
T = 11700°K , a = 0.1 (a ) , B = 0.0001 i (b„) . Solid lines: 
Asymptotic solution. Broken lines: Exact solution for same initial condi- 
tions. f. = t./(t, ], . 

1 1 '■ h'b 


FIG. 2. Inner solution. Argon, p^ = 1.5xl0~** atm., = (*^e^a ~ 

1000°K. T^ = 11700®K , = 0.1 * Solid lines. Asymptotic 

results for 8^ = 0.0001 i • Dotted lines: Asymptotic results for 

8 = 0.00059 (Sp) . Broken lines: Ref. 5 results for same initial 
conditions but with no excited level in atom model, u = u/Uj^ , T. = 

FIG. 3. Outer solution. Argon, p^ = 10~** atm., = 9610K, = 

5451°K , T^ = 19460°K , = 0.231 jf , B^ = 0.000219 i (Bj,)^ . 

Solid lines: Asymptotic solution. Dotted lines: Asymptotic solution for 
different initial conditions (Ref. 11). Broken lines: Exact solution for 
different initial conditions (Ref. 11). 


FIG. 4. Inner solution. Argon, p = lo”*^ atm., fT.l = 961®K, Tt 1 = 

a ' n' a '' e' a 

5451®K , T = 19460°K , a = 0.231 i (aj , B, = 0.000219 [b,,! . 

Solid lines: Asymptotic solution. Broken lines: Ref. 11 results for 
different atom model, p = p/p, , f. = T./fT,_l, . 

Jj 1 1 '• h'' b 

FIG. 5. Inner solution. Argon. Same initial conditions as in Fig. 4 but 

without heavy particle transport properties, p = p/p, , T, = T./Tt ) 

b i 1 ^ h"^ b 
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FIG. 6. Inner solution. Helium, p = 10 ^ atm., (t.) = 600®K , (t ) = 
4500®K , T = 23400°K , a = 0.0412 i (aj , 3 = 0.0000539 (bJ . 

C 3 3 3 3 

p = p/pj^ , T. = . 

FIG. 7. Inner solution. Argon, = 10 atm. , = 900°K , ~ 

6800°K , T = 22000°K , a = 0.205 ^ (aj , 3 = 0.000108 ^ (3p) . 

Solid lines: Asymptotic solution. Broken lines: Ref. 11 results for 

different atom model, p = p/p. , f. = T./fT, ), . 

D 1 1 '• n-'b 

FIG. 8. Inner solution. Argon. Same initial conditions as in Fig. 7 but 
without heavy particle transport properties, p = p/p-^ > ^i ~ '^i'^^'^h^b ‘ 

FIG. 9. Inner solution. Argon p^ = 10 ^ atm., ” 970®K , 

(t ) = 4600®K , T = 1700QOK , a = 0.250 / faJ , 3 = 0.000482 / 

e-'a ’ c ’a '• E-'a a 

(3p) . Solid lines: Asymptotic solution. Broken lines: Ref. 11 results 

for different atom model p = p/p-^ > ^i ” "^i^^^h^b ' 

FIG. 10. Inner solution. Argon. Same initial conditions as in Fig. 9 
but without heavy particle transport properties, p = p/p^ > 
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